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I.  INTRODUCTION 


A.  OVERVIEW 

Ships  can  be  initially  designed  for  a  variety  of  purposes  or  uses. 
But,  as  the  times  change,  they  may  be  no  longer  needed  for  this  original 
use  or  may  be  considered  for  other  types  of  operations  which  could 
require  it  to  operate  in  a  different  environment  than  originally  intended 
and  designed  for.  The  operating  characteristics,  in  this  new 
environment,  may  not  be  suitable  for  the  intended  use  of  the  ship. 

The  United  States  Air  Force  has  been  considering  the  use  of  T- 
AGOS  class  3  cmd  4  ships  as  sea-based  radar  platforms.  These  ships 
were  originally  designed  for  surveillance  and  ASW  operations.  The 
operating  characteristics  for  the  original  mission,  and  thus  what  the  ship 
was  designed  for,  are  different  to  what  would  be  required  for  this  current 
consideration.  To  act  as  a  sea-based  radar  platform,  consideration 
would  have  to  be  given  to  improving  the  seakeeping  qualities  of  the  ship, 
in  particular  with  regards  to  roll  motion.  Various  methods  exist  for 
improving  the  roll  of  ships  and  are  discussed  in  the  next  section. 

In  this  study,  the  effects  of  bilge  keels  on  the  ship  were  examined. 
In  particular,  various  sizes  of  bilge  keels  for  the  full  range  of  ship  speed 
and  sea  states.  A  FORTRAN  program,  SHIPMO  [Ref.  1],  was  utilized  to 
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perform  the  roll  angle  calculations  for  the  vaiying  conditions.  The  data 
obtained  was  utilized,  in  the  form  of  various  graphs  and  polar  plots,  to 
show  the  effects  of  bilge  keels.  Operabihty  indices  were  found  for  each  of 
the  bilge  keels  under  examination.  If  properly  compared,  these 
operability  indices  can  provide  a  useful  comparison  of  the  effectiveness  of 
bilge  keels. 


B.  STABILIZING  METHODS 


Severed  methods  exist  for  improving  the  seagoing  stability  of  the 
ship  and  include,  but  not  limited  to,  fin  stabilizers,  roll  tanks,  and  bilge 
keels  and  may  be  used  together  or  individually.  The  stability 
characteristics  of  the  ship  and  the  effects  of  employing  these  various 
methods  can  easily  be  seen  by  looking  at  the  roll  cheiracteristics  that  are 
experienced  by  the  ship.  Various  advantages  and  disadvantages  exist  for 
each  of  these  methods  and  must  be  considered  when  determining  a 
particular  method  to  be  used.  [Ref.  2] 

The  first  type  of  stabilization  technique  that  is  available  are  active 
fin  stabilizers.  These  produce  a  controlled  roll  moment  where  the  phase 
and  amplitude  is  such  that  it  counteracts  the  external  heel  moment.  The 
fins  are  most  effective  at  higher  speeds,  generally  greater  than  ten  knots, 
since  the  force  on  the  fin  varies  directly  proportional  to  the  speed  of  the 
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ship.  Below  ten  knots,  the  stabilizing  moment  available  is  not  adequate 
to  effectively  counter  the  heel  moment.  In  addition,  the  load  on  the  fin  or 
fin  pivot  may  be  such  that  the  fin  is  unable  to  undergo  full  angular 
displacement  at  certain  speeds.  This  could  result  in  degraded 
performance  at  other  additional  speeds.  Various  parameters  such  as 
shapes,  ship  locations,  and  angles  of  attack  can  be  varied  to  help  obtain 
the  maximum  performance  from  them.  The  following  are  advantages  and 
disadvantages  of  active  fin  stabilizers: 

Advantages: 

1.  They  offer  the  highest  possible  roll  reduction  with  no  reduction  in 
static  stability  characteristics.  They  are  the  most  effective  of  all  single 
stabilizing  devices. 

2.  They  are  used  in  ships  of  different  sizes. 

3.  They  inflict  very  small  increases  in  ship  resistance  and  have  small 
auxiMaiy  power  requirements. 

Disadvantages : 

1.  They  cire  not  effective  at  low  speeds. 
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2.  They  take  up  moderate  machinery  space,  especially  if  they  are 
retractable.  This  is  desirable  as  they  are  less  prone  to  damage. 

3.  High  initial  cost  due  to  the  controlling  equipment  and  machinery 
required. 

The  next  method  of  stabilization  that  exists  are  roll  tanks.  These 
tanks  can  be  active  or  passive  stabilizers  which  uses  ship’s  motion  in 
such  a  way  to  cause  water  in  the  tanks  to  move  in  a  direction  as  to 
oppose  the  ship’s  motion.  Two  different  tank  configurations  can  be 
employed  -  free  surface  or  U-tube.  A  major  limitation  to  this  method  is 
that  the  tanks  can  only  be  tuned  to  one  frequency.  This  is  usually  the 
natural  frequency  since  this  is  most  likely  where  the  largest  roll  angles 
will  occur.  They  are  not  as  effective  at  other  frequencies  and  can  actually 
increase  the  roll  angles  attained  at  these  other  frequencies.  More 
sophisticated  methods  exist  where  the  resistance  in  the  duct  can  be 
adjusted  for  the  frequency  of  the  exciting  waves  which  allows  dampening 
to  be  achieved  for  all  frequencies.  The  advantages  and  disadvantages  for 
this  method  of  roll  stabilization  are  as  follows; 

Advantages: 

1.  High  roll  reduction  rates  of  up  to  70%  are  possible. 
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2.  They  remain  effective  at  low  speeds. 


3.  Vary  little  auxiliaiy  power  is  required. 

4.  Moderate  initial  cost  and  low  maintenance  required. 

5.  Not  vulnerable  to  damage  since  they  are  internal  to  the  hull. 

Disadvantages: 

1.  Moderate  space  is  required  in  the  hull. 

2.  1-4%  reduction  in  deadweight  capacity  for  constant  displacement. 

3.  Reduction  in  initial  static  stability  due  to  free  surface  effects, 
especially  for  free  surface  tanks. 

The  fined  method  of  stabilization  that  was  considered  was  bilge 
keels.  Bilge  keels  are  the  most  widely  used  and  simplest  kind  of  roll 
stabilization  in  current  use.  They  consist  of  a  fin  fixed  to  the  hull  at  or 
near  the  bilge.  The  natural  period  of  the  roll  of  the  ship  is  proportional  to 
the  radius  of  gyration  of  the  ship.  By  attaching  this  fin,  the  radius  of 
gyration  of  the  ship  is  effectively  increased.  This  results  in  an  increased 
mass  of  water  to  roll  with  the  ship  and  therefore  an  increase  in  the 
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period  of  the  roll.  Under  forced  rolling  conditions,  such  as  in  a  seaway, 
the  increased  natural  period  that  results  from  the  bilge  keel  results  in  a 
roll  amplitude  reduction.  The  increased  resistance  to  roll  due  to  viscous- 
eddy  effects  supplied  by  the  bilge  keels  plays  an  even  larger  role  in  roll 
reduction.  Energy  is  dissipated  by  viscous  flow  from  around  the  ship 
and  this  energy  dissipation  is  increased  substantially  by  bilge  keel  use. 
The  advantages  and  disadvantages  of  this  method  are  as  follows: 

Advantages: 

1.  Bilge  keels  are  simple  eind  easy  to  fit. 

2.  They  remain  effective  at  relatively  low  speeds. 

3.  Negligible  reduction  in  ship’s  deadweight  capacity,  no  reduction  in 
initial  static  stability,  no  auxiliaiy  power  requirements,  and  negligible 
space  occupied  by  the  hull. 

4.  Low  initied  cost. 

Disadvantages : 

1.  Since  they  are  external  to  the  hull,  there  is  added  resistance  to  ahead 
motion  that  must  be  overcome  by  the  main  engines. 
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2.  Comparatively  to  other  methods,  bilge  keels  offer  smaller  amounts  of 
roll  reduction. 

3.  Bilge  keels,  carefully  aligned  to  flow  around  the  hull  in  calm  waters  to 
reduce  forward  motion  resistance,  can  lead  to  added  resistance  during 
roll  motion. 

4.  Vulnerable  to  damage. 

All  advantages  and  disadvantages  of  the  various  methods  must  be 
considered  when  determining  which  stabilization  method  is  best  suited 
for  the  application  of  interest.  Some  particular  advantages  or 
disadvantages  may  be  of  greater  importance  in  the  decision  process  than 
others.  Active  fin  stabilizers,  although  the  most  effective,  suffer  from 
severe  degradation  at  low  speeds.  Passive  roll  tanks  have  comparatively 
large  response  times  and  quick  changes  in  roll  motion  can  disrupt  the 
operation  of  the  roll  tanks.  Active  roll  tanks  require  large  initial  costs 
and  auxiliaiy  power  with  minimed  added  advantage.  Due  to  the 
simplicity  and  relative  effectiveness  at  low  speeds,  bilge  keels  were  the 
method  of  choice. 
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n.  EVALUATION  OF  RESPONSE 


A.  REGULAR  WAVE  MOTION 


An  important  consideration  for  any  floating  body  is  the  effects  of 
waves  on  the  free  surface,  especially  surface  waves  with  a  period  of  a  few 
seconds.  Other  waves,  such  as  subsurface  waves,  usually  have  much 
longer  periods  and  have  negligible  effect  compared  to  surface  waves. 
When  analyzing  these  types  of  waves,  some  assumptions  must  be  made. 
These  include  that  the  fluid  is  ideed  and  body  motions  are  sufficiently 
small  to  linearize.  From  this,  the  appropriate  fluid  mechanics  tools  can 
be  utilized  to  describe  sea  waves  and  ship  motions  based  on  potential 
flow  theory.  [Ref.  3] 

Potential,  or  ideal,  flow  theory  is  based  on  a  couple  of  assumptions. 
The  first  assumption  is  that  mass  is  preserved.  If  a  control  volume  is 
placed  around  the  object  of  interest  and  the  fluid  around  it,  then  the 
mass  which  enters  the  volume  must  either  accumulate  in  the  control 
volume  or  leave.  Utilizing  the  divergence  theorem,  conservation  of  mass, 
and  assuming  that  density  remains  constant,  the  final  form  of  the 
continuity  equation  inside  the  control  volume  becomes: 


du  ^dv  ^dw 
dx  dy  dz 


(1) 
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The  next  assumption  is  that  flow  is  irrotational.  This  results  in  the 


property  that  circulation  around  any  closed  curve  is  zero  or: 


Jt/Jr  =  0 


(2) 


From  equation  (2)  and  the  definition  of  a  velocity  potential,  the  continuity 
equation  can  be  written  as: 


vV  = 


d^<l> 


=  0 


(3) 


which  is  better  known  as  Laplace’s  Equation. 

Once  the  velocities  and  ^have  been  solved  for  and  assuming 
unsteady,  irrotational  flows,  pressure  can  be  found  with  the  help  of 
Bernoulli’s  equation  and  can  be  computed  from: 

/7  =  -p-^-^pV0-V^-pgz  (4) 

Tlie  plane  progressive  wave  system  is  the  simplest  free  surface 
wave  formation  scheme.  It  is  two  dimensional,  sinusoidal  in  time,  and 
propagates  with  a  phase  velocity  such  that  an  observer  moving  with  the 
same  velocity  will  make  the  wave  appear  stationary.  Adopting  a 
Cartesian  coordinate  system,  the  free  surface  elevation  cam  be  expressed 
in  the  following  form: 

rj[x,  t)  =  A  cos(kx  -  ox)  (5) 

where:  A-wave  amplitude 

k-wave  number 
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Ship  motions  induced  by  regular  wave  motion  are  of  most  concern. 
In  the  simplest  case,  the  waves  incident  upon  the  body  may  be  assumed 
as  plane  progressive  waves  of  small  amplitude  and  sinusoidal  time 
dependence.  Ship  motion  is  cdso  assumed  sufficiently  small  for  linear 
theory  to  hold.  Waves  incident  to  a  stationary  or  moving  body  will  cause 
the  body  to  move  with  six  degrees  of  freedom  or  to  heave,  pitch,  roll,  yaw, 
surge,  and  sway.  In  general,  the  body  motion  in  any  of  the  six  degrees  of 
freedom  can  be  expressed  as: 

=  (6) 

;=i 

where:  A- wave  amplitude 

F-exciting  force  per  unit  wave  amplitude 
D-matrix  of  equations  of  motion 

A  new  fundamental  quantity  can  also  be  introduced,  known  as  the 
Response  Amplitude  Operator  (RAO),  as  follows: 

Zj{Q),U,9)  =  ^  (7) 

A 

This  corresponds  to  the  complex  amplitude  of  body  motion  in  the  jth 
mode  in  response  to  an  incident  wave  of  unit  amplitude,  frequency,  and 
direction,  and  is  generally  known  as  the  transfer  function,  which  can  be 
calculated  once  the  added  mass,  damping,  exciting,  and  hydrostatic 
forces  are  known. 
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The  responses  derived  from  the  above  motions  are  important  for 
evaluating  the  seakeeping  performance  of  the  ship.  Computation  of  the 
RAO’s  is  simple  once  the  added  mass  and  damping  coefficients  are 
known.  But,  this  is  the  most  difficult  part.  A  technique  for  determining 
these  values  utilizes  the  Strip  theory.  In  this  method,  the  ship  is 
assumed  as  a  slender  body.  In  other  words,  the  beam  and  draft  are 
much  smaller  than  the  length  of  the  ship.  The  cross-section  of  the  ship 
are  divided  into  “strips”,  or  as  for  Frank’s  method  [Ref  4]-  a  series  of 
straight  line  segments  integrated  along  the  entire  length  of  the  ship  so  as 
to  determine  the  overall  effect.  By  utilizing  the  strip  theory  method,  the 
ship’s  motions  in  regular  waves  can  be  determined.  Through  a 
combination  of  the  above  equations,  the  resulting  linearized 
simultaneous  equations  that  must  be  solved  for  to  determine  these 
motions  are: 

+  +  +  J=1..6  (8) 

where:  M  -  mass  matrix 

A-added  mass  matrix 
B-dampenmg  matrix 
C-hydrostatic  restoring  force  matrix 

F'  -Froude  Krylov  exciting  force  in  the  jth  mode  of  motion 
F"  -DiSraction  exciting  force  in  the  jth  mode  of  motion 
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-complex  amplitude  of  motion  in  the  kth  direction 

B.  ROLL  DAMPING 

Roll  damping  can  be  accomplished  by  a  variety  of  methods  as 
discussed  above.  The  roll  damping  hydrodynamic  moment  used  in 
equations  of  motion  tends  to  be  nonlinear  in  nature  with  respect  to  the 
roll  angle.  To  allow  roll  damping  to  be  calculated,  as  in  SHIPMO  [Ref  5], 
it  must  first  be  linearized.  Utilizing  the  method  described  by  Himeno 
(1981)  [Ref  6],  the  roU  damping  coefficient  can  be  assumed  in  regular 
waves  as: 

^44(^4  )  =  -Si  +  4^2  (9) 

5k 

where:  B  -  damping  matrix/coefficients 

C  -  Complex  amplitude  of  motion 
ft)  -  wave  frequency 


The  damping  coefficients  are  a  function  of  frequency,  the  point 
about  which  roll  occurs,  and  the  forward  speed.  In  regular  wave  motion 
with  constant  forward  speed,  the  damping  coefficient  values  are  constant 
about  a  fixed  point.  When  irregular  or  random  waves  are  encountered, 
the  following  equation  by  Himeno  applies: 


where: 


^44(^4)  -  A  +  j~^^2 

V  TT 

O’  -  variance  in  angular  roll  velocity 


(10) 
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The  damping  coefficients  are  composed  of  the  following 
components: 


where: 


^2  ~  ^BK 


(11) 


Bp  -  frictional  damping  component 

-  wave  generation  damping  component 

-  hull  lift  damping  component 

B  ^  -  hull  eddy  damping  component 


Bpf.  -  bilge  keel  damping  component 
The  above  components  are  based  mainly  on  experimental  analysis 
which  are  measured  about  a  fixed  point  (roll  center)  and  are  usually 
limited  in  their  range  of  applicability.  If  the  beam  to  draft  ratios  are 
much  more  than  about  2.5,  the  eddy  dampening  by  Himeno  tends  to 
overpredict  the  eddy  component.  In  addition,  determining  the  point 
about  which  the  ship  is  rolling,  or  roll  center,  can  prove  to  be  difficult.  In 
particular,  a  ship  in  a  seaway  has  only  an  instantaneous  roU  center  and 
can  be  calculated  as  follows: 


RC  = 


If 


|cos(a2  -0C4)-  tan{(oj  +  0:4)  sin(a2  -  0^4 )} 


(12) 


where: 


0^2  -  sway  phase  angle 


-  roll  phase  angle 
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C.  MOTIONS  IN  A  SEAWAY 


Wave  patterns  in  an  open  sea  are  ever  changing  with  respect  to 
time  and  space,  in  a  manner  that  appears  to  defy  analysis  whether  it  is 
linear  or  second  order  Stokes.  [Ref  7  ]  Ambient  waves  on  the  surface  of 
the  sea  are  dispersive  as  well  as  random  in  nature  where  random  refers 
to  the  character  of  the  wave  height  distribution.  The  distribution  of 
sinusoidal  waves  are  such  that  they  have  continuously  distributed 
amplitudes  and  phases  such  that  the  summation  of  the  variation  of  wave 
height  with  time  is  not  systematic  in  any  respect,  but  random.  The 
generating  mechanism  is,  predominantly,  the  effect  of  wind  in  the 
atmosphere  upon  the  water  surface.  Spectral  Density,  S(©),  provides 
useful  data  that  has  been  removed  from  a  random  wave  record,  h(t).  The 
random  wave  record  can  be  processed  such  that  a  S(®)  vs.  wave 
frequency,  g),  curve  can  be  generated.  The  spectral  density  is  obtained 
from  a  wave  height  record  taken  during  a  period  of  time  when  sea 
conditions  are  considered  unchanging  (stationaiy)  for  a  certain  sea  state. 
The  spectral  energy  density,  S{(o,e),  (or  directional  energy  spectrum)  can 
be  integrated  over  all  wave  directions,  6 ,  to  give  the  frequency  spectrum 
as  given  below: 

2n 

S{(0)  =  ^S{(o,d)dd  (13) 
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For  open  seas,  it  is  appropriate  to  assume  that  waves  are  unidirectional  which 
allows  the  energy  spectrum  to  be  proportional  to  a  delta  function  in  0  where  the  wave 
crests  can  be  considered  parallel  and  the  fluid  motion  to  be  two  dimensional.  Wave 
spectra  of  this  form  is  considered  long  crested  and  is  sufficient  to  describe  the  wave 
enviroranent. 

In  naval  architecture  and  ocean  engineering,  it  is  usually  appropriate  to  assume 
waves  to  be  long  crested.  With  this  simplification,  existing  information  for  the  energy 
spectrum,  based  on  theory  and  full  scale  observations,  can  be  used  for  the  energy 
spectrum.  Usually,  we  are  more  concerned  with  larger  waves.  The  most  common 
parameter  that  takes  this  into  account  is  the  significant  wave  height,  Hy^ ,  which  is 
defined  as  the  average  of  the  highest  one  third  of  all  waves  and  can  be  expressed  as 
follows: 

;/,„=4.0K)”  (14) 

Where  is  the  integral  of  S(co )  over  all  frequencies.  More  generally,  the  moments, 

nif ,  of  the  spectrum  can  defined  by: 

00 

rrij  =  Ja)'iS(a))Jto,/ =  0,1,2,...  (15) 

0 

Also,  the  average  frequency  of  the  spectrum,  defined  as  the  expected  number  of  zero 
upcrossings  of  the  body  that  is  considered  per  unit  time  (i.e.  the  number  of  times  the 
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wave  amplitude  passes  through  zero  with  a  positive  slope),  can  be 
determined  as  follows: 


m,  =  (^)‘” 


(16) 


The  average  period  between  zero  upcrosstngs  is: 


(17) 


Figure  1 .  Typical  Pierson-Moskowltz  wave  spectrum 
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A  good  model  for  fuUy  developed  seas  is  the  Pierson-Moskowitz 
spectrum.  Figure  1  shows  a  typical  Pierson-Moskowitz  spectrum.  This 
spectrum  is  based  on  significant  wave  height  where  it  is  assumed  that 
the  wave  spectra  has  reached  a  steady-state  equilibrium  and  is 
independent  of  duration  and  fetch.  In  other  words,  the  waves  are 
considered  fully  developed.  It  is  recognized  as  an  asymptotic  form 
reached  only  sifter  an  extended  period  of  wind  with  no  contamination 
from  an  underl3dng  swell  and  can  be  given  as  follows: 

8(0))  =  0.008  Ig ^  /  6)^  exp|-0.032(g  /  H^O)^  j  (18) 

where:  g  -  acceleration  of  gravity 

-  significant  wave  height 
CO  -  wave  frequency 

Once  the  frequency  spectrum  is  found,  another  quantity  can  be 
determined  which  will  allow  the  spectral  description  of  waves  to  be 
generalized  into  regular  harmonic  waves.  As  previously  discussed,  the 
Response  Amplitude  Operator,  RAO,  is  a  nondimensional  quantity  that 
can  be  used  to  relate  the  amount  of  heave  amplitude  per  unit  wave 
length.  If  the  sea  waves  are  described  by  a  rsindom  distribution  and  if 
body  response  to  each  wave  is  defined  by  a  RAO,  the  body  response  will 
be: 
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where: 


11,(0  =  9?JjzX®,e>'“’'^i4(co,e) 


(19) 


91  -  any  body  response 
Z  -  response  amplitude  operator 
CO  -  wave  frequency 
9  -  wave  direction 

The  body  response  is  also  a  random  variable.  The  seaway  spectrum  may  be 
related  to  any  body  response  by  the  following  equation: 

‘S'«(<o)=K(®)|25(ffl)  (20) 

where:  '  spectrum  of  the  seaway 

The  above  RAO  is  not  only  valid  for  regular  wave  motions  but  also  in  a 
spectrum  of  random  waves.  This  allows  for  the  assumption  that  a  vessel  with  favorable 
characteristics  in  regular  waves  will  also  have  the  same  favorable  characteristics  in 
irregular  waves.  From  equation  20,  various  ship’s  motions  can  be  calculated  from  a 
given  sea  spectrum,  such  as  the  Pierson-Moskowitz  discussed  previously. 
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D.  CALCULATION  PROGRAM 

The  computer  program  SHIPMO  [Ref  1]  was  utilized  to  make  the 
various  calculations,  more  specifically  roll  calculations,  needed  for  this 
study.  The  computer  program  calculates  a  variety  of  motion  responses 
and  is  able  to  predict  ship  motions  in  six  degrees  of  freedom.  Regular 
wave  motions  are  calculated  utilizing  a  modified  strip  theory  of  Salvesen, 
Tuck,  and  Faltinsen  (1970)  [Ref  8]  which  is  based  on  a  slender  body  in 
an  ideal  fluid  with  corrections  made  for  viscous  and  other  effects  .  Two 
dimensional  properties  are  calculated  using  Frank’s  method  (1967)  [Ref 
9]  with  a  lid  applied  to  remove  any  irregulEU*  frequencies  which  can  cause 
calculation  difficulties  if  not  applied.  Frank’s  method  is  computationally 
fast  and  any  ship’s  cross-section  can  be  approximated  with  as  much 
accuracy  as  desired.  The  wave  spectrum  can  either  be  user  inputted  or 
chosen  from  one  of  six  methods  available  in  the  program.  All 
calculations,  for  this  study,  were  based  on  the  classical  Pierson- 
Moskowitz  spectrum  as  described  in  the  previous  section. 

The  program  reads  in  the  initial  ship  data  and  calculates  the 
potential  for  the  four  modes  of  motion  -  pitch,  roll,  yaw,  and  sway,  over 
the  specified  frequency  range  at  each  designated  cross  section  utilizing 
Frank’s  method.  The  calculations  are  based  at  the  ships’  amidships 
centerline  position  unless  otherwise  specified.  The  progreun  allows  for  a 
maximum  of  21  stations  to  be  selected  along  the  length  of  the  ship,  with 
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a  maximum  of  15  input  points  for  each  station.  These  are  required  to  provide  the  hull 
shape  data  necessary  for  the  various  calculations.  The  wave  spectral  ordinates  are 
calculated  by  the  Pierson-Moskowitz  method  as  discussed  above.  The  spectral 
ordmates  are  also  used  to  determine  the  response  amplitude  operator  as  was  also 
previously  discussed.  The  motion  amplitude  is  calculated  in  all  six  degrees  of  freedom 
with  the  horizontal  and  vertical  plane  motions  solved  separately.  Roll  dampening  is 
calculated  by  the  Himeno  method  with  eddy  dampening  for  box  barges  with  sharp 
comers  being  calculated  by  either  Ikeda’s  (1977)  [Ref  10]  or  Yamashita  and  Katagiri’s 
method  (1980).  [Ref  11]  For  bilge  keel  damping,  the  Himeno  method  is  also  utilized. 
From  this,  the  desired  ships  motions  can  be  calculated  over  ship’s  speed,  heading 
angle,  and  wave  frequency  variations. 

The  program  required  a  data  initialization  file,  SHIPMO.IN,  and  is  included  in 
Appendix  A.  The  file  inputs  a  variety  of  constants  required  for  the  proper 
computations  to  be  made.  Bilge  keel  width,  depth,  and  location  are  inputted  from  this 
file.  Station  positions  make  up  the  majority  of  the  file  which  include  the  station 
position  relative  to  the  center  of  the  ship.  For  each  station,  hull  depth  and  width 
positions  relative  to  the  center  and  water  lines,  are  provided  with  the  first  input  point 
located  at  the  keel  of  the  ship. 
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E.  SHIP  AND  ENVIRONMENTAL  CONSIDERATIONS 


The  analysis  performed  in  this  study  were  based  on  the  Class  4  T-AGOS  ship. 
Table  1  provides  some  of  the  general  characteristics  of  the  ship  that  were  necessary  for 
the  calculations. 


Length 

224  feet 

Location  of  Centerline 

97  feet  aft  of  bow 

Design  Draft 

15  feet 

Maximum  Ship  Speed 

11.3  knots 

Table  1.  General  Characteristics  T-AGOS  Class  4  Ships 


SHIPMO,  as  described  in  the  previous  section,  was  utilized  to  determine  the 
roll  characteristics  of  the  ship  based  on  different  bilge  keel  characteristics.  Hull  width 
and  depths  for  various  locations  along  the  length  of  the  ship  were  required  for  the 
program’s  calculations.  These  were  obtained  from  the  body  plans  for  the  T-AGOS  4 
class  ship.  [Ref  12]  A  sealed  version  of  the  one  actually  used  may  be  viewed  in 
Appendix  B.  The  actual  data  inputted  for  each  of  the  reference  points  may  be  found  in 
the  data  input  file  SHIPMO. IN  in  Appendix  A. 

Roll  calculations  were  made  through  the  entire  speed  and  sea  state  spectrum. 
Speeds  of  the  ship  were  incremented  from  zero  to  12  knots  in  two  knot  speed 
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increments.  Sea  states  are  based  on  and  provide  an  indicator  of  the 
heights  of  the  waves,  usually  the  significant  wave  height,  being 
encountered  by  a  vessel.  They  range  firom  1,  relatively  calm  seas,  to  8 
which  are  usually  only  encountered  during  strong  storms  such  as 
hurricanes.  Table  2  provides  the  significant  wave  height  for  each  of  the 
sea  states  considered. 


Sea  State 

Significant  Wave  Height  (feet) 

1 

1.9 

2 

4.1 

3 

5.7 

4 

7.4 

5 

13.0 

6 

20.8 

7 

40.3 

8 

61.6 

Table  2.  Sea  State  vs.  Significant  Wave  Height 
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m.  ROLL  RESPONSE  RESULTS 


A.  ROLL  CALCULATIONS 

As  discussed,  SHIPMO  was  utilized  to  make  the  roll  calculations 
necessary  for  this  study  and  a  sample  input  file  is  included  in  Appendix 
A.  The  ship’s  roll  characteristics  were  analyzed  based  on  several  bilge 
keel  sizes  including  no  bilge  keel.  Appendix  C  provides  a  sample  of  the 
program  that  was  utilized  to  analyze  the  data  obtained  and  print  out  the 
results.  Table  3  lists  the  different  bilge  keel  sizes  that  were  ainalyzed. 
The  maximum  bUge  keel  width  was  based  upon  two-thirds  the  length  of 
the  ship  and  subsequently  divided  in  half  for  each  smaller  size.  A 
maximum  width  of  three  feet  was  set  and  reduced  in  size  in  one  foot 
increments.  Each  bilge  keel  was  analyzed  such  that  it  was  assumed 
centered  along  the  centerline  of  the  ship. 
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To  ensure  that  a  good  representative  data  set  was  obtained  for  each  bilge  keel 
configuration,  enough  data  points  were  required.  For  each  ship’s  speed  analyzed,  the 
roll  angle  was  found  for  each  sea  state  at  22.5  degree  increments  of  wave  direction 
around  a  360  degree  circle  centered  on  the  ship. 

Once  this  data  was  obtained,  a  curve  of  significant  roll  angles  placed  on  the  ship 
for  each  sea  state  and  a  given  speed  and  bilge  keel  configuration,  could  be  found 
utilizing  the  spline  curve  fitting  Ofiinction  of  MATLAB  [Ref  13]  for  each  degree  of 
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wave  motion  from  0  to  360  degrees.  Enough  data  was  required  to 
ensure  a  smooth  fit  of  the  curves.  To  illustrate  this  point,  a  graph  of  the 
roll  angles  for  no  bilge  keel  and  a  ship’s  speed  of  2  knots  may  be  seen  in 
Figure  2.  Similar  graphs  for  all  the  other  ship’s  speeds  and  bilge  keel 
configurations  could  also  be  produced  if  needed. 


Roll  angle  vs.  wave  angle  (2  kts.,  sea  states  (ss)  1  to  8) 


ss8 

ss7 

ss6 

ss5 

ss4 

ss3 

ss2 

ssl 


Figure  2.  Roll  Angle  vs.  wave  angle  (no  bilge  keel  and  2  knots  ship  speed) 
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From  the  data  obtained  as  just  discussed,  a  spltned  curve  of  roll 
angles  for  all  sea  states  for  each  degree  of  wave  motion  was  found.  From 
this,  the  maximum  allowed  sea  states  to  produce  a  roll  angle  of  2,  4,  or  6 
degrees  were  found.  A  polar  plot  of  this  data  for  each  bilge  keel 
configuration  and  respective  speeds  are  shown  in  figures  3  to  72.  In 
addition,  polar  plots  combining  all  ship  speeds  for  a  certain  bilge  keel 
configuration  were  generated  and  can  be  observed  in  figures  73  to  102. 
The  next  section  discusses  the  effectiveness  of  bilge  keels.  Operability 
indices  were  also  determined  and  plotted  to  help  show  this. 

By  examining  the  figures,  several  observations  become  evident.  As 
was  expected,  the  maximum  cdlowed  sea  state  increased  as  the 
maximum  allowed  roll  angle  was  increased  and  as  ship’s  speed  was 
increased,  the  maximum  allowed  sea  state  correspondingly  decreased. 
Beam  seas  seemed  to  have  the  least  effect  on  the  roll  characteristics  of 
the  ship  and  roll  was  assumed  zero  for  sea  directions  relating  to  the  bow 
or  stem  of  the  ship.  The  largest  effect  on  roll  emgle  came  from  quartering 
seas  with  quartering  seas  aft  the  beam  having  the  largest  effect.  Thus, 
aft  quartering  seas  provided  the  limiting  position  for  roll  angle  limits.  As 
can  be  seen,  significant  roll  reduction  can  be  expected  by  the  use  of  bilge 
keels. 
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Figure  3.  Sea  state  vs.  wave  angle  for  no  bilge  keel  -  0  kts. 


Figure  4,  Sea  state  vs,  wave  angle  for  no  bilge  keel  -  2  kts. 


Maximum  sea  state  vs.  wave  angle  (for  8  kt.  and  2,4,6  degree  roll  angle) 
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Figure  7.  Sea  state  vs.  wave  angle  for  no  bilge  keel  -  8  kts. 


Maximum  sea  state  vs.  wave  angle  (for  10  kt.  and  2,4,6  degree  roll  angle) 
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Maximum  sea  state  vs.  wave  angle  (for  2  kt.  and  2,4,6  degree  roll  angle) 
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Figure  11.  Sea  state  vs.  wave  angle  for  150x3  bilge  keel  -  2  kts 

Maximum  sea  state  vs.  wave  angle  (for  4  kt.  and  2,4,6  degree  roll  angle) 
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Figure  12.  Sea  state  vs.  wave  angle  for  150x3  bilge  keel  -  4  kts. 


Maximum  sea  state  vs.  wave  angle  (for  6  kt.  and  2,4,6  degree  roll  angle) 
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Figure  13.  Sea  state  vs.  wave  angle  for  150x3  bilge  keel  -  6  kts. 


Figure  15.  Sea  state  vs.  wave  angle  for  150x3  bilge  keel  -  10  kts 


Maximum  sea  state  vs.  wave  angle  (for  12  kt.  and  2,4,6  degree  roll  angle) 
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Figure  16.  Sea  state  vs.  wave  angle  for  150x3  bilge  keel  -  12  kts. 


Maximum  sea  state  vs.  wave  angle  (for  0  kt.  and  2,4,6  degree  roll  angle) 
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Figure  17.  Sea  state  vs.  wave  angle  for  150x2  bilge  keel  -  0  kts. 


Maximum  sea  state  vs.  wave  angle  (for  2  kt.  and  2,4,6  degree  roll  angle) 
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Figure  18.  Sea  state  vs.  wave  angle  for  150x2  bilge  keel  -  2  kts. 


Figure  19.  Sea  state  vs.  wave  angle  for  150x2  bilge  keel  -  4  kts. 


Maximum  sea  state  vs.  wave  angle  (for  6  kt.  and  2,4,6  degree  roll  angle) 
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Figure  20.  Sea  state  vs.  wave  angle  for  150x2  bilge  keel  -  6  kts. 
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Figure  21.  Sea  state  vs.  wave  angle  for  150x2  bilge  keel  -  8  kts 
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Figure  22.  Sea  state  vs.  wave  angle  for  150x2  bilge  keel  -  10  kts. 


Figure  23.  Sea  state  vs.  wave  angle  for  150x2  bilge  keel  -  12  kts. 


Figure  24.  Sea  state  vs.  wave  angle  for  150x1  bilge  keel  -  0  kts. 


Figure  25.  Sea  state  vs.  wave  angle  for  150x1  bilge  keel  -  2  kts. 
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Figure  26.  Sea  state  vs.  wave  angle  for  150x1  bilge  keel  -  4  kts. 


Figure  27.  Sea  state  vs.  wave  angle  for  150x1  bilge  keel  -  6  kts 
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Figure  28.  Sea  state  vs.  wave  angle  for  150x1  bilge  keel  -  8  kts. 


Figure  29.  Sea  state  vs.  wave  angle  for  150x1  bilge  keel  -  10  kts. 
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Figure  30.  Sea  state  vs.  wave  angle  for  150x1  bilge  keel  -  12  kts. 


Maximum  sea  slate  vs.  wave  angle  (for  0  kt.  and  2,4,6  degree  roll  angle) 

90„ 


-  2  degrees 
*  4  degrees 
‘  6  degrees 


Figure  31.  Sea  state  vs.  wave  angle  for  75x3  bilge  keel  -  0  kts 
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Figure  32.  Sea  state  vs.  wave  angle  for  75x3  bilge  keel  -  2  kts. 
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Figure  33.  Sea  state  vs.  wave  angle  for  75x3  bilge  keel  -  4  kts. 


Figure  34.  Sea  state  vs.  wave  angle  for  75x3  bilge  keel  -  6  kts 


Figure  35.  Sea  state  vs.  wave  angle  for  75x3  bilge  keel  -  8  kts. 


Figure  36.  Sea  state  vs.  wave  angle  for  75x3  bilge  keel  -  10  kts. 


Figure  37.  Sea  state  vs.  wave  angle  for  75x3  bilge  keel  -  12  kts 


Maximum  sea  state  vs.  wave  angle  (for  0  kt.  and  2,4,6  degree  roll  angle) 


Figure  38.  Sea  state  vs.  wave  angle  for  75x2  bilge  keel  -  0  kts. 


Figure  41.  Sea  state  vs.  wave  angle  for  75x2  bilge  keel  -  6  kts 


Figure  42.  Sea  state  vs.  wave  angle  for  75x2  bilge  keel  -  8  kts 


Figure  43.  Sea  state  vs.  wave  angle  for  75x2  bilge  keel  -  10  kts. 


Maximum  sea  state  vs.  wave  angle  (for  1 2  kt.  and  2,4,6  degree  roll  angle) 


Figure  44.  Sea  state  vs.  wave  angle  for  75x2  bilge  keel  -  12  kts. 


Figure  45.  Sea  state  vs.  wave  angle  for  75x1  bilge  keel  -  0  kts. 


Figure  46.  Sea  state  vs.  wave  angle  for  75x1  bilge  keel  -  2  kts. 


Figure  47.  Sea  state  vs.  wave  angle  for  75x1  bilge  keel  -  4  kts 


Maximum  sea  state  vs.  wave  angle  (for  6  kt.  and  2,4,6  degree  roll  angle) 


270 

Figure  48.  Sea  state  vs.  wave  angle  for  75x1  bilge  keel  -  6  kts 


Figure  49.  Sea  state  vs.  wave  angle  for  75x1  bilge  keel  -  8  kts 


Maximum  sea  state  vs.  wave  angle  (for  10  kt.  and  2,4,6  degree  roll  angle) 
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Figure  50.  Sea  state  vs.  wave  angle  for  75x1  bilge  keel  -  10  kts 
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Figure  51.  Sea  state  vs.  wave  angle  for  75x1  bilge  keel  -  12  kts. 
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Figure  52.  Sea  state  vs.  wave  angle  for  37.5x3  bilge  keel  >  0  kts. 
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Figure  53.  Sea  state  vs.  wave  angle  for  37.5x3  bilge  keel  -  2  kts 
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Figure  54.  Sea  state  vs.  wave  angle  for  37.5x3  bilge  keel  -  4  kts 
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Figure  55.  Sea  state  vs.  wave  angle  for  37.5x3  bilge  keel  -  6  kts. 
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Figure  56.  Sea  state  vs.  wave  angle  for  37.5x3  bilge  keel  -  8  kts. 
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Maximum  sea  state  vs.  wave  angle  (for  0  kt.  and  2,4,6  degree  roll  angle) 

90„ 


270 

Figure  59.  Sea  state  vs.  wave  angle  for  37.5x2  bilge  keel  -  0  kts. 
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Figure  61.  Sea  state  vs.  wave  angle  for  37.5x2  bilge  keel  -  4  kts. 
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Figure  62.  Sea  state  vs.  wave  angle  for  37.5x2  bilge  keel  -  6  kts. 
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Maximum  sea  state  vs.  wave  angle  (for  8  kt.  and  2,4,6  degree  roll  angle) 
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Figure  63.  Sea  state  vs.  wave  angle  for  37.5x2  bilge  keel  -  8  kts. 
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Figure  64.  Sea  state  vs.  wave  angle  for  37.5x2  bilge  keel  -  10  kts. 
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Maximum  sea  state  vs.  wave  angle  (for  2  kt.  and  2,4,6  degree  roll  angle) 
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Figure  67.  Sea  state  vs.  wave  angle  for  37.5x1  bilge  keel  -  2  kts. 


Figure  69.  Sea  state  vs.  wave  angle  for  37.5x1  bilge  keel  -  6  kts. 


Figure  70.  Sea  state  vs.  wave  angle  for  37.5x1  bilge  keel  -  8  kts. 


Figure  71.  Sea  state  vs.  wave  angle  for  37.5x1  bilge  keel  -  10  kts. 
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Figure  72.  Sea  state  vs.  wave  angle  for  37.5x1  bilge  keel  -  12  kts. 
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B.  OPERABILITY  INDICES 


Operability  Indices,  as  used  in  this  study,  can  provide  a  relative 
measure  of  the  ship’s  ability  to  operate  under  all  sea  states  and  speeds. 
From  the  polar  plot  for  each  speed  and  a  given  maximum  roll  angle,  the 
area  encircled  by  the  resulting  curve  could  be  calculated  utilizing  the 
trapezoidal  method  by  the  use  of  the  trapz  function  in  MATIAB. 

Similarly,  the  maximum  area  of  the  circle  can  also  be  calculated  if  the 
ship  is  assumed  to  be  able  to  operate  in  all  sea  states  without  exceeding 
the  maximum  allowed  roll  angle.  The  operability  index  can  then  be 
calculated  by  dividing  the  calculated  area  of  interest  by  the  maximum 
area.  Appendix  D  provides  the  program  used  to  calculate  the  operability 
indices. 

By  comparing  the  operability  index  related  to  a  particular  bilge  keel 
size  and  ship’s  speed  to  that  of  no  bilge  keel,  a  relative  measure  of  the 
effectiveness  of  that  particular  bilge  keel  can  be  determined.  Figure  103 
is  a  plot  of  the  Operability  indices  for  various  ship’s  speeds,  no  bilge  keel, 
and  2,4,  and  6  degree  maximum  roll  angles.  Since  one  objective  of  this 
study  was  to  determine  the  effectiveness  of  various  bilge  keels  on  the  roll 
characteristics  of  the  ship,  a  better  measure  of  the  effectiveness  was 
found  by  presenting  all  bilge  keel  operability  indices  with  respect  to  the 
no  bilge  keel  operability  indices.  In  other  words,  the  operability  indices 
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of  the  bilge  keel  under  consideration  were  divided  by  the  respective 
operability  indices  for  the  ship  with  no  bilge  keel.  If  the  two  values  were 
exactly  the  same,  the  result  would  be  1.0.  By  multiplying  the  result  with 
the  corresponding  no  bilge  keel  operability  index  value,  one  could 
determine  the  actual  respective  operability  index  value.  Figures  104  to 
112  are  plots  of  these  relative  values  for  the  range  of  ship’s  speed. 

Further  examination  of  the  operability  indices  provided  some 
useful  points.  First,  the  operability  indices  seem  lowest  at  2  and  8  knots 
for  the  no  bilge  keel  ship  with  the  absolute  minimum  at  8  knots.  Six 
knots  and  speeds  above  ten  knots  provide  the  speeds  best  suited  to  allow 
the  ship  the  highest  amount  of  operability.  This  results  in  the  bilge  keels 
having  less  of  an  effect  at  these  speeds  eind  the  greatest  effect  at  2  and  8 
knots.  As  can  be  seen,  any  bilge  keel  arrangement  still  significantly 
improves  the  operability  of  the  ship  at  most  any  ship  speed.  It  appears 
that  shortening  the  length  of  the  bilge  keel  by  half  while  maintaining  the 
same  width  reduces  the  effectiveness  by  about  0.2.  If  the  length  is 
maintained  constant  while  the  width  is  reduced  in  size  by  one  foot,  the 
overall  reduction  in  effectiveness  appears  to  be  approximately  0.2  also. 
These  observations  are  generalities  and  can  vaiy  by  a  greater  or  less 
extent  based  on  the  different  bilge  keels  and  roll  angles  examined. 
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Operability  Index  for  no  bilge  keel  and  2,4,  and  6  deg  roll  angle 
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Figure  103.  Operability  Index  for  no  bilge  keel. 


Operability  Index  for  150x3  bilge  keel  and  2,4,  and  6  deg  roll  angle 
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Figure  104.  Operability  Index  for  150x3  bilge  keel. 
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Operability  Index  Comparison  Operability  Index  Comparison 


Operabilify  Index  Comparison  Operability  Index 


Figure  107.  Operability  Index  for  75x3  bilge  keel. 


Figure  108.  Operability  Index  for  75x2  bilge  keel. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 


As  was  seen,  not  only  the  size  of  a  wave  incident  to  a  ship  but  also 
the  ship’s  speed  can  produce  large  roU  effects.  Wave  direction  also 
played  a  large  part  on  the  magnitude  of  the  ship’s  roll.  The  data  provided 
favorable  results  regarding  the  use  of  bilge  keels  to  reduce  the 
magnitude  of  the  ship’s  roll.  These  results,  though,  are  based  on  theory 
and  historical  results.  So,  the  possibility  exists  that  actual  bilge  keel 
performance  could  and  probably  will  vary  from  the  results  found  in  this 
study.  The  following  conclusions  can  be  drawn: 

1.  As  was  expected,  the  maximum  allowed  sea  state  increased  as  the 
maximum  allowed  roll  angle  was  increased. 

2.  As  ship’s  speed  was  increased,  the  maximum  allowed  sea  state 
correspondingly  decreased. 

3.  Beam  seas  seemed  to  have  a  significant  effect  on  the  roll 
characteristics  of  the  ship  and  roll  was  zero  for  sea  directions  relating  to 
the  bow  or  stem  of  the  ship.  The  largest  effect,  however,  on  roll  angle 
came  from  quartering  seas  with  quartering  seas  aft  the  beam  having  the 
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largest  effect.  Thus,  ail  quartering  seas  being  the  limiting  position  for 
roll  angle  limits. 

4.  Significant  roll  dampening  can  be  applied  by  the  use  of  bilge  keels. 

5.  The  operability  indices  seem  lowest  at  2  and  8  knots  for  the  no  bilge 
keel  ship  with  the  absolute  minimum  at  8  knots.  Thus,  bilge  keels  will 
have  the  greatest  effect  at  these  speeds. 

6.  Any  bilge  keel  arrangement  can  significantly  improve  the  operability  of 
the  ship. 

7.  Speeds  below  six  knots  and  above  ten  knots  provided  the  ship  the 
highest  amount  of  operability  with  no  bilge  keel.  Bilge  keels  still  proved 
to  significantly  decrease  the  roll  motion  of  the  ship  even  at  these  speeds. 

8.  Shortening  the  length  of  the  bilge  keel  by  half  while  maintaining  the 
same  -width  appeared  to  reduce  the  effectiveness  by  approximately  0.2. 

9.  If  the  bilge  keel  length  is  maintained  constant  while  the  width  is 
reduced  in  size  by  one  foot,  the  overall  reduction  in  effectiveness  appears 
to  be  approximately  0.2  also. 
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B.  RECOMMENDATIONS 


To  continue  the  research  on  this  topic,  the  following 
recommendations  are  made: 

1.  This  study  investigated  the  effects  of  bilge  keels  on  the  roll 
characteristics  of  the  ship.  Several  other  methods  of  providing  roll 
dampenmg  exist  and  could  also  be  explored. 

2.  Once  the  placement  of  the  radar  masts,  or  other  equipment  of 
concern,  have  been  determined,  motion  characteristics  at  these  locations 
could  be  determined.  These  motions  would  result  not  only  because  of 
roll  but  also  because  of  the  other  modes  of  motion  of  the  ship. 

3.  Once  the  maximum  allowable  roll  angle  for  proper  operation  of  the 
masts  and/ or  equipment  of  concern  has  been  determined,  the  minirmim 
required  bilge  keel  size  can  then  be  determined  by  the  same  methods  of 
this  study. 
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APPENDIX  A.  SAMPLE  SHIPMO  DATA  INPUT  FILE. 
Initial  for  T-agos 

0012101100100  17  0 

224.00  30.00  3.00  3.00  22.96  5.28  12.00 

224.0000  1.9905  32.1740  O.OOOOOE+00  O.OOOOOE+00  0.0000 

18.7500  -18.7500  1.0000 

1  97.0000  0.0000  0 

0.0000  -15.0000 
6  87.3000  0.0000  0 

0.0000  -15.0000 
0.5208  -12.0000 
1.8750  -9.0000 
3.0458  -6.0000 
5.0000  -3.0000 
6.5000  0.0000 

6  77.6000  0.0000  0 

0.0000  -15.0000 
3.1563  -12.0000 
5.8125  -9.0000 

8.5417  -6.0000 

9.9375  -3.0000 
10.9790  0.0000 

6  67.9000  0.0000  0 

0.0000  -15.0000 
5.4479  -12.0000 
10.0000  -9.0000 
12.4896  -6.0000 
13.6250  -3.0000 
14.7083  0.0000 

6  58.2000  0.0000  0 

0.0000  -15.0000 
7.9792  -12.0000 
14.4375  -9.0000 
15.6146  -6.0000 
16.6875  -3.0000 

17.5417  0.0000 

6  48.5000  0.0000  0 

0.0000  -15.0000 
11.0729  -12.0000 
17.0000  -9.0000 
17.9479  -6.0000 
19.1667  -3.0000 
19.5625  0.0000 

6  38.8000  0.0000  0 
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0.0000  -15.0000 
14.4896  -12.0000 

18.6250  -9.0000 
19.5313  -6.0000 

20.6250  -3.0000 

20.8750  0.0000 

6  19.4000  0.0000  0 

0.0000  -15.0000 
19.0938  -12.0000 
19.9375  -9.0000 
20.6354  -6.0000 

21.5000  -3.0000 

21.5000  0.0000 

6  0.0000  0.0000  0 
0.0000  -15.0000 
19.5833  -12.0000 
20.0625  -9.0000 
20.8229  -6.0000 

21.5000  -3.0000 

21.5000  0.0000 
6  -19.4000  0.0000  0 

19.5833  -15.0000 
20.1875  -12.0000 
20.8229  -9.0000 

21.5000  -6.0000 

21.5000  -3.0000 

21.5000  0.0000 
6  -38.8000  0.0000  0 

0.0000  -15.0000 
19.1354  -12.0000 

20.1250  -9.0000 

20.7500  -6.0000 

21.5000  -3.0000 

21.5000  0.0000 

6  -58.2000  0.0000  0 

0.0000  -11.6250 

19.1250  -9.0000 

19.6250  -7.5000 
20.1146  -6.0000 

20.5000  -3.0000 

21.5000  0.0000 

6  -67.9000  0.0000  0 

0.0000  -9.2500 

3.1250  -9.0000 

13.7500  -7.5000 
19.1042  -6.0000 


94 


19.9375  -3.0000 
21.4792  0.0000 

6  -77.6000  0.0000  0 

0.0000  -6.8750 
8.1667  -6.0000 

18.2500  -4.5000 

19.2500  -3.0000 

20.3750  -1.5000 

21.3750  0.0000 

6  -87.3000  0.0000  0 

0.0000  -4.7500 
3.0000  -4.6250 
6.0000  -4.2500 

17.2500  -3.0000 

18.8750  -1.5000 
20.6250  0.0000 

6  -97.0000  0.0000  0 

0.0000  -3.0000 
6.0000  -2.7500 
12.0000  -2.5000 
15.0000  -2.3750 
18.0000  -0.6250 
18.9792  0.0000 

1  -103.8759  0.0000  0 

0.0000  -2.0000 
5.2800  0.0000 

9999.0000  0.0000  0.000 

1.0000  0.2000  2.0000  0.1000  13.5025  13.5025  0.0000 

0.0000  180.0000  22.5000 
1.9000 
0.0 
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APPENDIX  B.  BODY  PLANS  FOR  T-AGOS  CLASS  4  SHIP 
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APPENDIX  C.  SAMPLE  DATA  AND  OUTPUT  FILE. 


%  Sets  up  required  data  sets  for  sea  state,  roll  angle,  and  angle  on  ship 

%  Speed  =  0-12  kt.,  Sea  State  =  1:8.  150x1  bilge  keel 

clear 

thl=[0.0  22.5  45.0  67.5  90.0  112.5  135.0  157.5  180.0]; 
ss=[l  2  3  4  5  6  7  8]; 
speed=[0  2  4  6  8  10  12]; 

w=[0.0000,0. 1 190,0.2590,0.3350,0. 1930,0.2490,0.2340,0. 1390,0.0000; 

0.0000,0.8550,1.4100,1.3900,0.8300,1.1400,1.0700,0.7080,0.0000; 

0.0000,1.5700,2.3600,2.2400,1.5200,1.9400,1.9100,1.2900,0.0000; 

0.0000,2.2900,3.2400,3.0500,2.3100,2.7500,2.7700,1.9600,0.0000; 

0.0000,3.7500,4.9600,4.7500,4.0600,4.5000,4.5600,3.4200,0.0000; 

0.0000,4.5800,6.0000,5.8200,5.1800,5.5700,5.5900,4.2600,0.0000; 

0.0000,5.2600,6.8700,6.7600,6.2800,6.6100,6.5100,4.9200,0.0000; 

0.0000,5.4400,7.2000,7.1700,6.7300,7.0300,6.8500,5.1400,0.0000; 

0.0000,0.1760,0.3400,0.3730,0.1900.0.2220,0.1790,0.0991,0.0000; 

0.0000,1.1700,1.7200,1.5000,0.8200,1.0300,0.8610,0.5190,0.0000; 

0.0000,2.0000,2.7400,2.3700,1.5000,1.7600,1.5700,0.9640,0.0000; 

0.0000,2.7400,3.6000,3. 1900,2.2800,2.5500,2.3800, 1 .5600,0.0000; 

0.0000,4.1000,5.2000,4.8100,4.0400,4.3400,4.2300,3.0800,0.0000; 

0.0000,4.8200,6.1200,5.8200,5.1300,5.5000,5.3900,4.0000,0.0000; 

0.0000,5.3400,6.9400,6.7000,6.2300,6.5400,6.4100,4.7800,0.0000; 

0.0000,5.5300,7.2400,7.0800,6.6800,7.0300,6.8100,5.0400,0.0000; 

0.0000,0.2900,0.4660,0.4180,0.1860,0.2000,0.1420,0.0740,0.0000; 

0.0000,1.6500,2.1000,1.6000,0.7990,0.9320,0.7070,0.3970,0.0000; 

0.0000,2.5600,3. 1500,2.4800, 1.4400, 1.5800, 1.3000,0.7390,0.0000; 

0.0000,3.2700,3.9900,3.2500,2.1900,2.3100,2.0100,1.2300,0.0000; 

0.0000,4.4400,5.4700,4.8100,3.9200,4.1000,3.8400,2.6900,0.0000; 

0.0000,5.0300,6.3000,5.7600,5.0600,5.3100,5.0800,3.7100,0.0000; 

0.0000,5.4800,6.9700,6.6000,6.1100,6.4800,6.2200,4.5600,0.0000; 

0.0000,5.6700,7.2300,6.9600,6.5700,6.9600,6.6600,4.8600,0.0000; 

0.0000,0.5900,0.7010,0.4760,0.1840,0.1810,0.1150,0.0572,0.0000; 

0.0000,2.3800,2.6000,1.7200,0.7860,0.8600,0.6050,0.3190,0.0000; 

0.0000,3.2100,4.5700,2.5500, 1.3700, 1.4400, 1. 1000,0.5940,0.0000; 

0.0000,3.8100,4.3000,3.2500,2.0300,2.0800,1.7000,0.9630,0.0000; 

0.0000,4.7500,5.6300,4.6600,3.6400,3.7700,3.3800,2.2300,0.0000; 

0.0000,5.2400,6.3400,5.5300,4.7700,5.0400,4.7000,3.3000,0.0000; 

0.0000,5.6200,7.0000,6.4000,5.8800,6.2600,5.9300,4.2500,0.0000; 

0.0000,5.7500,7.2300,6.7500,6.3600,6.7900,6.4200,4.6100,0.0000; 

0.0000,1.5500,1.2400,0.5460,0.1810,0.1650,0.0959,0.0455,0.0000; 

0.0000,3.5700,3.4500,1.9100,0.7850,0.8000,0.5270,0.2640,0.0000; 

0.0000,4.2700,4.3500,2.7400,1.3600,1.3500,0.9770,0.5000,0.0000; 

0.0000,4.6900,4.9800,3.4000, 1 .9900, 1 .9600, 1 .5200,0.8 180,0.0000; 
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0.0000,5.4300,6.0700,4.6800,3.5200,3.5700,3.1200,1.9700,0.0000; 
0.0000,5.7700,6.7000,5.4700,4.6100,4.8200,5.6200,3.0200,0.0000; 
0.0000,6.0700,7.2300,6.2400,5.7300,6.1400,5.7400,4.0100,0.0000; 
0.0000,6.1700,7.4400,6.5700,6.2100,6.6700,6.2800,4.3900,0.0000; 
0.0000,0.9200,2.5200,0.6370,0. 1 780,0. 1 500,0.08 1 1 ,0.0372,0.0000; 
0.0000,3.1600,4.7300,2.1400,0.7810,0.7440,0.4640,0.2230,0.0000; 
0.0000,3.9800,5.4900,2.9900,1.3700,1.2700,0.8690,0.4280,0.0000; 
0.0000,4.4900,6.0100,3.6600,2.0000,1.8700,1.3700,0.7070,0.0000; 
0.0000,5.2800,6.9400,4.8600,3.5300,3.4900,2.9500,1.8100,0.0000; 
0.0000,5.6700,7.4400,5.5700,4.6000,4.7700,4.2800,2.8700,0.0000; 
0.0000,5.9700,7.8600,6.2800,5.7100,6.1200,5.6500,3.9200,0.0000; 
0.0000,6.0700,7.9800,6.6000,6.2100,6.6700,6.2300,4.3300,0.0000; 
0.0000,0.5980,2.3000,0.7630,0.1750,0.1380,0.0698,0.0310,0.0000; 
0.0000,2.0400,5.2400,2.4300,0.7770,0.6950,0.4140,0.1930,0.0000; 
0.0000,2.7200,6. 1000,3.2900, 1.3600, 1 . 1900,0.7790,0.3720,0.0000; 
0.0000,3.2100,6.6900,3.9200,2.0000,1.7700,1.2400,0.6170,0.0000; 
0.0000,4.0400,7.5300,5.0600,3.5400,3.4000,2.7800,1.6600,0.0000; 
0.0000,4.4900,8.0000,5.7700,4.6100,4.7200,4.1700,2.7500,0.0000; 
0.0000,4.8400,8.3900,6.4400,5.7000,6.1200,5.6000,3.8600,0.0000; 
0.0000,4.9600,8.5500,6.7500,6.2000,6.6900,6.2100,4.2900,0.0000]; 


%  Utilizes  a  spline  function  to  smooth  data 

forj=  1:7 

ctr=l; 

en=(j-l)*8  +  1; 
beg=en+7; 
for  i=beg:-l:en 
t=0:180; 

wl  (i,  :)=:spline(th  1  ,w(i, :)  ,t); 
wlp(ctr,:)=wl(i,:): 
ctr=ctr+l; 
end 

figure(j) 

plot(t,wlp(;,:)) 

xlabelCWave  Angle  (degrees  -  150x1  bilge  keel)') 

ylabelCRoll  Angle  (degrees)') 

s=speed(i); 

titled'Roll  angle  vs.  wave  angle  (',int2str(s),'  kts.,  sea  states  (ss)  1  to 

8)']) 

Iegend('ss8','ss7','ss6','ss5',’ss4','ss3','ss2','ssl',-1) 

end 
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%  Calculates  maximum  roll  angle  (2,4,6  degrees)  for  angle  on  ship. 

ctrl=8; 

ctr2=l; 

for  j  1=1:7 

beg=(il-l)*8  +  1; 

en=beg+7; 

ctr=l; 

for  k=beg:en 

wl  a(ctr,  :)=w  1  (k, :) ; 
ctr=ctr+l; 

end 

ra2=8.0*ones(l,181); 
ra4=8 .  O*ones(  1,181); 
ra6=8 .  O*ones(  1 , 1 8 1) ; 

[rl,cl]=size(wla); 
ssn=  1 .00:0.0 1 :8.00; 
ssnc=size(ssn,2): 

for  i=l:cl 

wn=splme(ss,w  1  a(:  ,i) ,  ssn) ; 
if  wn(l,l)  >  2.0 
ra2(i)=0.0; 

end 

if  wn(l,l)  >  4.0 
ra4(i)=0.0; 

end 

if  wn(l,l)  >  6.0 
ra6(i)=0.0; 

end 

for  j=l:(ssnc-l) 

if  wn(l,j)  <=  2.0 
k=j+l; 

if  wn(l,k)  >  2.0 

rei2(i)=ssn(j); 

end 

end 

if  wn(l,j)  <=  4.0 
k=j+l; 

if  wn(l,k)  >  4.0 

ra4(i)=ssn(j): 

end 
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end 


if  wn(l,j)  <=  6.0 

k=j+l; 

if  wn(l,k)  >  6.0 

ra6(i)=ssn(j); 

end 

end 

end 

end 


%  Sets  up  roll  angle  and  theta  matrices  for  plotting 

fori=l:181 

z=i+180; 

th3(i)=t(i)*pi/180; 

th3(z)=t(i)*pi/180  +  pi; 

j=182-i; 

ra2(z)=ra2(j); 

ra4(z)=ra4(j); 

ra6(z)=ra6(j); 

end 


%  Prints  out  polar  plot  for  sea  state  limits  for  roll  angles  of  2,4,  and 

%  6  degrees  for  each  speed. 

fLgure(ctrl) 

polar(th3 ,  ra6 ,  'r') 

hold  on 

polar(th3  ,ra4,  'b') 
polar(th3 ,  ra2 , 'g') 

legend('g','2  degrees','b',’4  degrees', 'r','6  degrees',- 1) 

xlabelC  150x1  bilge  keel') 

hold  off 

s2  =  speed(jl); 

titled'Maximum  sea  state  vs.  wave  angle  (for  ',int2str(s2),'  kt.  and  2,4,6 

degree  roll  angle)']) 

pra2(ctr2,  :)=ra2(  1 , :); 

pra4(ctr2 ,  :)=ra4(  1 , :) ; 

pra6(ctr2,  :)=ra6(  1 ,:); 

ctrl=ctrl+l; 

ctr2=ctr2+l; 

end 

%Prints  out  polar  plots  for  sea  state  limits  for  aU  speeds  for  2,4,  and 
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%6  degrees  roll  angles. 
figure(15) 

polar(th3,pra2(  1 , ,  'y') 
legendCO  kts’,-1) 
hold  on 

polar(th3  ,pra2(2, :) , 'm') 
polar(th3,pra2(3, :) ,  'c') 
polar(th3,pra2{4, :) ,  'r') 
polar(th3,pra2(5, :) ,  'g') 
polar(th3,pra2(6, :), 'b') 
polar(th3,pra2(7, :) ,  'w') 
xlabelC  150x1  bilge  keel') 

titleCMaximum  sea  state  vs.  wave  angle  (for  all  speeds  and  2  degree  roll 
angle)') 

legendCy’.'O  kts',’m','2  kts',’c’,'4  kts','b’,'10  kts’,'w','12  kts',-1) 
hold  off 


fLgure(16) 

polar(th3,pra4(  1 , :) ,  'y') 
hold  on 

polar(th3,pra4{2, :) ,  m') 
polar{th3,pra4(3, :), 'c') 
polar(th3,pra4(4, :) ,  'r') 
polar(th3,pra4(5, :) , 'g') 
polar(th3  ,pra4(6, :) ,  'b') 
polar(th3,pra4(7, :),  'w') 
xlabelC  150x1  bilge  keel') 


titleCMaximum  sea  state  vs.  wave  angle  (for  all  speeds  and  4  degree  roll 
angle)') 

legend('y',’0  kts','m','2  kts','c','4  kts','b','10  kts',’w','12  kts',-1) 
hold  off 


figure!  17) 

polar(th3,pra6(  1  ,:),'y') 
hold  on 

polar(th3  ,pra6(2, :) , 'm') 
polar(th3,pra6(3, :), 'c') 
polar(th3,pra6(4, :) ,  'r') 
polar(th3,pra6(5,:),'g') 
polar(th3,pra6(6, :) ,  'b') 
polar(th3,pra6(7, :)  ,'w') 
hold  off 

xlabelC  150x1  bilge  keel') 

titleCMaximum  sea  state  vs.  wave  angle  (for  all  speeds  and  6  degree  roll 
angle)') 

legend('y’,'0  kts','m','2  kts','c','4  kts','b','10  kts','w','12  kts',-1) 
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figure(l) 

print  -depsc2  gl50xl0.eps 
figure(2) 

print  -depsc2  gl50xl2.eps 
figure(3) 

print  -depsc2  gl50xl4.eps 
figure(4) 

print  -depsc2  gl50xl6.eps 
fignre(5) 

print  -depsc2  gl5Qxl8.eps 
figure(6) 

print  -depsc2  gl5Qxll0.eps 
figiire(7) 

print  -depsc2  gl50xll2.eps 
figure(8) 

print  -depsc2  pl50xl0.eps 
figure(9) 

print  -depsc2  pl50xl2.eps 
figure(lO) 

print  -depsc2  pl50xl4.eps 
figure(ll) 

print  -depsc2  pl50xl6.eps 
figiire(12) 

print  -depsc2  pl50xl8.eps 
figured  3) 

print  -depsc2  pl50xll0.eps 
figure(14) 

print  -depsc2  pl50xll2.eps 
figure(15) 

print  -depsc2  pl50xl2d.eps 
figure(16) 

print  -depsc2  pl50xl4d.eps 
figure(17) 

print  -depsc2  pl50xl6d.eps 

%  Calculates  areas  under  2,4,  and  6  degree  curves  to  find  operability 

index 

for  i=l:7 

a  1 50 1 2(i)=trapz(th3  ,pra2(i, :)) ; 
a  1 50 1 4(i)=trapz(th3  ,pra4(i, :)) ; 
a  1 50 1 6(i)=trapz(th3,pra6(i, :)) ; 

end 

diary  on 
al5012 
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al5014 
al5016 
diary  off 


APPENDIX  D.  OUTPUT  PILE  FOR  OPERABILITY  INDICES 


%  Operability  index  for  no  and  all  bilge  keel  sizes 


sp=lO  2  4  6  8  10  12]; 

oi=[20.5352  19.6437  20.9513  21.8156  18.7591  22.2139 


23.1490; 

27.5689 

30.2916 

40.6990 

50.2655 

50.2655 

31.5754 

50.1601 

50.2655 

23.6220 

33.8953 

43.4954 

33.1574 

50.2655 

50.2655 

27.0449 

41.0418 

50.2655 

22.0652 

30.9436 

36.6973 

28.3462 

43.5711 

50.2655 

24.0074 

34.6294 

45.0302 

21.3010 

29.3208 

33.6042 


25.4916 

28.4353 

40.6389 

50.2655 

50.2655 

31.4809 

50.1964 

50.2655 

23.6932 

33.9282 

43.5624 

33.1825 

50.2655 

50.2655 

26.8404 

40.9971 

50.2655 

22.1496 

30.9119 

36.9689 

28.0363 

43.6137 

50.2655 

24.0688 

34.7334 

45.0246 

21.2152 

28.9641 

33.5519 


26.3032 

29.7059 

40.4951 

50.2655 

50.2655 

31.0794 

49.9981 

50.2655 

23.5707 

34.0430 

43.7010 

33.0098 

50.2655 

50.2655 

26.5161 

40.9294 

50.2655 

22.1688 

30.7803 

37.3480 

27.6600 

43.4905 

50.2655 

23.9180 

34.8465 

45.2023 

21.3695 

28.7246 

33.9582 


27.0467 

30.6955 

40.1558 

50.2655 

50.2655 

30.4644 

49.5953 

50.2655 

23.6384 

33.5065 

44.2738 

32.4087 

50.2655 

50.2655 

26.0630 

41.1273 

50.2655 

22.0700 

30.6469 

38.3693 

27.1524 

43.4775 

50.2655 

23.6660 

35.1848 

45.7350 

21.3614 

28.6705 

34.9195 


24.5390 

28.4974 

39.2916 

50.2655 

50.2655 

29.0074 

49.2371 

50.2655 

22.6258 

32.3755 

43.9998 

30.8072 

50.2655 

50.2655 

26.1045 

41.0348 

50.2655 

21.2421 

29.1149 

37.8269 

26.0138 

42.7571 

50.2655 

22.7242 

33.7725 

45.5974 

20.2741 

27.1758 

33.6370 


27.0156 

30.8023 

36.7968 

50.2655 

50.2655 

28.0754 

47.6217 

50.2655 

22.2212 

32.3211 

43.0580 

29.5798 

49.3792 

50.2655 

24.8517 

38.2164 

50.2655 

21.4560 

29.6535 

37.2617 

25.6165 

41.5179 

50.2655 

22.4861 

33.0551 

44.6895 

21.2836 

28.4457 

33.9735 


28.5117; 

32.1364; 

35.2270; 

50.2655; 

50.2655; 

28.4607; 

46.2013; 

50.2655; 

24.0342; 

33.3536; 

42.8258; 

29.5969; 

47.5442; 

50.2655; 

25.7132; 

38.5906; 

49.1335; 

23.5047; 

31.1143; 

37.9979; 

26.4812; 

40.6086; 

50.1873; 

24.2154; 

33.9994; 

44.3698; 

23.1232; 

30.0060; 

35.0874]; 


%  changes  operability  index  to  %  for  no  bilge  keel 
max=50.2665; 


for  i=l:3 

pooi(i,  :)=oi(i, :) .  /max; 
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end 


%  compares  operability  index  of  no  bilge  keel  to  different  bilge  keels 

for  i=4:3:30 
l=i+ 1 
m=i+2; 
for  j=l:7 

pooi(i,j)=oi(i  j) .  /oi(  1  ,j); 
pooi(l,j)=ol(l,j).  /oi(2,j); 
pooi(m  J)=oi(m,j) .  /oi(3,j)  ; 

end 

end 


%Printing  out  results 
flgure(l) 

plot(sp,pooi(l,:),sp,pooi(2,:),sp,pooi(3,:)) 
xlabelf'Speed  (kts)’) 
ylabelCOperability  Index') 

titleCOperability  Index  for  no  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legend{'2  deg','4  deg','6  deg',-1) 
print  -depsc2  oinbk.eps 
figure(2) 

plot(sp,pooi(4,:),sp,pooi(5,:),sp,pooi(6,:)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  comparison') 

titleCOperability  Index  for  150x3  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legendC2  deg','4  deg','6  deg',-1) 
print  -depsc2  oil50x3.eps 
figure(3) 

plot(sp,pooi(7,:),sp,pooi(8,:),sp,pooi(9,:)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  Comparison') 

titleCOperability  Index  for  150x2  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legend('2  deg','4  deg','6  deg',-1) 
print  -depsc2  oil50x2.eps 
figiu:e(4) 

plot(sp,pooi(10,:),sp,pooi(ll,:),sp,pooi(12,:)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  Comparison') 

titleCOperability  Index  for  150x1  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legend('2  deg','4  deg','6  deg',-1) 
print  -depsc2  oil50xl.eps 
figure(5) 
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plot(sp,pooi(13,:),sp,pooi(14,:),sp,pooi{15,:)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index') 

titleCOperabilily  Index  for  75x3  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
Iegend('2  deg', '4  deg', '6  deg',-1) 
print  -depsc2  oi75x3.eps 
figure(6) 

pIot(sp,pooi(16,:),sp,pooi(17,:),sp,pooi(18,:)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  Comparison') 

titleCOperability  Index  for  75x2  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
IegendC2  deg', '4  deg', '6  deg',-1) 
print  -depsc2  oi75x2.eps 
figure{7) 

pIot(sp,pooi(  19, :)  ,sp,pooi(20, :)  ,sp,pooi{2 1 , :)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  Comparison') 

titleCOperability  Index  for  75x1  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legend('2  deg','4  deg','6  deg',-1) 
print  -depsc2  oi75xl.eps 
fLgure(8) 

plot(sp,pooi(22, :)  ,sp,pooi(23, :)  ,sp,pooi(24, :)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  Comparison') 

titleCOperability  Index  for  37.5x3  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legendC2  deg','4  deg','6  deg',-1) 
print  -depsc2  oi375x3.eps 
figure(9) 

plot(sp,pooi(25, :)  ,sp,pooi(26, :)  ,sp,pooi(27, :)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  Comparison') 

titleCOperability  Index  for  37.5x2  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legendC2  deg','4  deg','6  deg',-1) 
print  -depsc2  oi375x2.eps 
figure!  10) 

plot(sp,pooi(28, :) ,  sp  ,pooi(29, : )  ,sp,pooi(30, :)) 
xlabelCSpeed  (kts)') 
ylabelCOperability  Index  Comparison') 

titleCOperability  Index  for  37.5x1  bilge  keel  and  2,4,  and  6  deg  roll  angle') 
legend('2  deg','4  deg','6  deg',-1) 
print  -depsc2  oi375xl.eps 
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